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Abstract. Recent radioastronomical observations of Fara- 
day rotation in the solar corona can be interpreted as evi- 
dence for coronal currents, with values as large as 2.5 x 10 9 
Amperes (ISpanglerL 120071) . These estimates of currents are 
used to develop a model for Joule heating in the corona. It 
is assumed that the currents are concentrated in thin current 
sheets, as suggested by theories of two dimensional magne- 
tohydrodynamic turbulence. The Spitzer result for the resis- 
tivity is adopted as a lower limit to the true resistivity. The 
calculated volumetric heating ra te is compared with a n inde- 
pendent theoretical estimate by ICranmer et al ( 2007 ). This 
latter estimate accounts for the dynamic and thermodynamic 
properties of the corona at a heliocentric distance of several 
solar radii. Our calculated Joule heating rate is less than the 
Cranmer et al estimate by at least a fact or of 3 x 10 5 . The cur- 



rents inferred from the observations of Spa ngler I (12007h are 
not relevant to coronal heating unless the true resistivity is 
enormously increased relative to the Spitzer value. However, 
the same model for turbulent current sheets used to calculate 
the heating rate also gives an electron drift speed which can 
be comparable to the electron thermal speed, and larger than 
the ion acoustic speed. It is therefore possible that the coro- 
nal current sheets are unstable to current-driven instabilities 
which produce high levels of waves, enhance the resistivity 
and thus the heating rate. 



1 Introduction 



In a recent paper, Spangler (2007) reported radioastronomi- 
cal observations which were consistent with the presence of 
coronal currents in the range of hundreds of MegaAmperes 
to a few GigaAmperes. This measurement was made using 
Faraday rotation observations of a radio source occulted by 



the corona, and the coronal plasma probed was at heliocentric 
distances of 5.2 to 6.7 R Q . In the present paper, I discuss the 
implications of these observations for the process of coronal 
heating by Joule heating. 



As discussed in I Spangler (2007) the currents reported 
(and summarized in Section 2 below) correspond to the net 
current within an Amperian loop defined by the two, closely- 
spaced lines of sight through the corona to the different parts 
of the radio source. The measured net current could be, and 
probably is, a residual due to numerous current filaments 
with alternate positive and negative current density within 
the Amperian loop. 

This topic is of interest because Joule heating has been 
identified as the primary mec hanism for heating the c losed- 



field part of the corona dGudiksen and Nordlund , 2005 



Peter et all 120061) . The purpose of this paper is to make 
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model-dependent estimates of the heating rate due to Joule 
dissipation of these currents. As expected, the calculation 
involves introduction of several "imponderables", i.e. phys- 
ical characteristics of the turbulence in the corona which are 
poorly constrained by observations, but which play an im- 
portant role in coronal heating. I feel this exercise is worth- 
while in identifying coronal parameters which are important 
in coronal heating, so that they can be targeted for future ob- 
servational investigations. 

The outline of this paper is as follows. In Section 2, 
I briefl y summarize the observational results of ISpanglerl 
(2007) which were the basis of the estimates of the coronal 
current. Section 3 is the most important part of the paper; it 
introduces a model for current-carrying coronal turbulence, 
and identifies the most important characteristics of this tur- 
bulence. A glossary of the variables and parameters intro- 
duced in this discussion is given in Table 1. This model is 
used to obtain an estimate of the volumetric heating rate due 
to Joule heating. Section 4 briefly considers the possibility 
that current densities in these sheets could be large enough to 
generate turbulence via current-driven instabilities, and thus 
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produce high wave levels which enhance the resistivity and 
thus the Joule heating rate. Section 5 summarizes what has 
been learned from this exercise and presents conclusions. 



2 Brief Summary of Radioastronomical Measurements 
of Coronal Currents 



The result reported by ISpangler I ((2007) was of a difference 
in the Faraday rotation measure ARM between two lines 
of sight to two components of an extragalactic radio source 
(3C228). These lines of sight were separated by an angu- 
lar distance 9, which corresponds to a linear separation in 
the corona between the two lines of sight, I = 9d, where 
d is the distance to the Sun. In the observations reported 

46 arcseconds and I = 33, 000 



Spangler (120071) 9 



m 

km. The observations were made when the line of sight to 
3C228 passed through the corona at heliocentric distances 
from 5.2 — 6.7-R f^. The technique is illustrated in Figure 1 of 
Spangler I d2007h . 

The fundamental physical relation used in t he technique is 
expressed by equation (3) of ISpangler I d2007l) 



ARM = C <h nB ■ ds ~ Cn <h B ■ ds 



(1) 



where the integral is around an Amperian loop through the 
corona, consisting of the two lines of sight, closed by imag- 
inary line segments which join the two lines of sight, at lo- 
cations infinitely separated from the corona. In this formula, 
C is a collection of atomic constants which arise in the de- 
scription of Faraday rotation, defined as C = 87r ^ c Ti eom 2 = 
2.631 x 10~ 13 in SI units, n(x) is the electron density, B(x) 
is the vector magnetic field in the corona, and ds is an incre- 
mental step around the Amperian Loop. 
Use of Ampere's Law 



B ■ ds = /jqI 



(2) 



in equation (1) shows that the differential rotation measure 
ARM is directly related to the current within the Amperian 
loop defined by the two lines of sight. The transition from 
the middle to the right term in equation (1) involves an ap- 
proximation, in which the position-dependent plasma den- 
sity in the integrand is replaced by an effective me an density 
n. This approximation is discussed at length in I Spangler I 



(2007), where arguments for its plausibility are presented. 

As is the case with Ampere's Law, the differential Fara- 
day rotation measurement provides information on the net 
current within the Amperian Loop. In general, we expect 
both positive and negative currents to be flowing within the 
loop. The electrical current flowing in a given location in 
the corona could be much larger than that deduced by the ar- 



Spangler I (120071) presented the following results from two 
observing sessions with the Very Large ArrajQ radiotele- 
scope. Each session lasted approximately 8 hours. 

1. In one of the two sessions, there was a confident de- 
tection of a ARM event, with a corresponding inferred 
electrical current of 2.5 x 10 9 Amperes. 

2. In the second observing session, a marginal ARM 
event was detected with a value of / (which may well 
be considered an upper limit) of 2.3 x 10 8 Amperes. 

3. During a several hour period of good data quality, no 
significant ARM events were detected, with a corre- 
sponding upper limit to the current of 7.7 x 10 8 Am- 
peres. 

4. Although the data from th e earlier investigation of 



Sakurai and Spangler I (119941) have not been reanalysed 



in this manner, examination of Figure 1 1 of that paper 
shows no clear evidence of a ARM event in several 
more hours of VLA observation. All of this indicates 
that detection of a clear differential rotation measure 
event between lines of sight separated by ~ 30, 000 km 
is relatively rare. 

The aforementioned observations will represent the obser- 
vational constraints imposed on the theory developed in the 
next section. 



3 Implications for Coronal Heating 

In this section, I discuss the implications of the results from 
Spangler I (120071) for coronal heating. The presence of elec- 
trical currents indicates that Joule heating will occur as well. 
I will calculate an estimate of the average volumetric heating 
rate of a system of currents which could produce the obser- 
vations discussed in Section 2. 

This calculation will be highly model dependent, as well 
as dependent on assumptions regarding the nature of the 
current sheets. Since the subsequent discussion will intro- 
duce many assumed parameters of the coronal current sheets, 
some of which are poorly constrained, I will follow the once- 
common practice in physics and astronomy literature of in- 
cluding a glossary of physical variables. This is contained 
in Table 1 . The following analysis assumes that the current 
is contained in a number of thin current sheets within the 
Amperian Loop. This model is illustrated in Figure 1. A 
coordinate system is defined by having one axis (the z axis) 
coincide with that of the large scale coronal magnetic field. 
I assume that the current sheets are extended along the large 
scale field, as is the case in quasi-2D magnetohydrodynamics 



gumen ts above, and contained in equation (7) of ISpangler 
(2003). 



The Very Large Array is an instrument of the National Ra- 
dio Astronomy Observatory. The NRAO is a facility of the Na- 
tional Science Foundation, operated under cooperative agreement 
with Associated Universities, Inc. 
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Variable 


Definition 


A 


Area of the Amperian Loop defined by the experiment 




ratio of current sheet width to domain size, L c /A 


£ 


Joule heating per current sheet (Watts) 


E 


Joule heating per domain, due to all current sheets 


e 


Domain-averaged volumetric heating rate 


£T 


Domain-averaged volumetric heating rate for turbulent current sheets 


es 


Domain-averaged volumetric heating rate for deterministic current sheets 


n 


Electrical resistivity of a plasma, — 1/a 


r]s 


Spitzer resistivity (equation (16)) 


f+ 


Fraction of sheets with positive current density within Amperian Loop 


f- 


Fraction of sheets with negative current density within Amperian Loop 


I 


Electrical current per sheet 


hbs 


Total current within the Amperian Loop 


j 


Current density within one current sheet 


I 


Transverse separation of two lines of sight in the corona 


L c 


Width of a current sheet 


A 


Transverse scale of a domain which contains a few current sheets 


A z 


Extent of domain along the large scale magnetic field 


/' 


Anisotropy of the current sheets (scale length along magnetic field/scale length across field) 


N 


Number of current sheets within one domain 


N+ 


Number of sheets with positive current density within Amperian loop 


N- 


Number of sheets with negative current density within Amperian loop 


N T 


Total number of current sheets within the Amperian loop defined by the experiment 


Slos 


Effective line-of-sight depth of the coronal plasma 


a 


electrical conductivity of a plasma 


to 


Thickness of a current sheet 


vr> 


electron drift speed 


V 


Volume of a current sheet 


Vn 


Volume of a domain 


Z c 


Extent of a current sheet along the large scale magnetic field 



dZank and Matfhaeus 1 119921) . The current sheet properties 
have a weaker dependence on the coordinate along the large 
scale field than on the coordinates in a plane perpendicular 
to that field. I begin by assuming that one can define a "do- 
main" which has a scale A perpendicular to the large scale 
coronal magnetic field, and which contains a small integer 
number N of current sheets. In the analysis which follows in 
Sections 3.1 and 3.2, 1 will assume that all current sheets are 
identical. The current sheet properties which are introduced 
are obviously to be understood as mean values from a distri- 
bution. Figure 1 illustrates such a domain. The extent of the 
domain in the direction perpendicular to the plane defined by 
Figure 1 is A 2 = /iA, with fi > 1. The current sheets have 
a width L c , a. thickness t c , and an extension along the large 
scale field Z c > L c . 

The picture which has been draw n so far is consistent with 
the original view of lParker ( 1972 ). It is also consistent with 
results from studies of 2D MHD turbulence, which show 
that turbulent evolution results in the formation of isolated, 
intense sheets of current and vorticity. The deve l opmen t 
which follows is based on results from I Spangler I ((1999), 



which contains an extensive bibliography to the literature 
where these ideas were dev eloped earlier, most importantly 



Zank and Matthaeus (1992). 



In the case of current sheets which arise from 2D MHD 
turbulence, the number of positive and negative current 
sheets should be equal, and the expectation value of the cur- 
rent in an Amperian loop is zero. The detection of net cur- 
rents (via differential Faraday rotation) would then be inter- 
preted as a statistical fluctuation of the total current about the 
zero expectation value. In what follows, I will refer to these 
as "turbulent current sheets". 

It is also possible that the physics of the corona selects 
current sheets with a preferred sign of the current density, at 
least for that portion of the corona which is probed in a Fara- 
day rotation experiment. This situation is referred to as that 
of "deterministic current sheets", and is discussed in Section 
3.2. Within this model, I assume that the properties of the 
individual current sheets are essentially the same as in the 
turbulent model, but that there is a preference for one sign 
of current density. It should be noted that the turbulent cur- 
rent sheet model is based on analytic and numerical solutions 
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Fig. 1. A vision of current-carrying coronal turbulence. The 2 di- 
mensional plane represented is perpendicular to the radial direction, 
which is also the direction of the large scale coronal magnetic field. 
The figure portrays the current density in gray-scale format. Black 
represents large positive current density, white is large negative cur- 
rent density, and gray indicates zero current density. This diagram 
illustrates the basic model for the current in the corona. It is con- 
tained within intense, narrow current sheets of both positive and 
negative sign. There are a few such current sheets, with thickness 
t c and width L c (noted in the fi gure as "t" and "L "), within a "do- 
main" of width A. Adapted from Spangler] dl999l) . 



of the equations of 2D magnetohydrodynamics, whereas the 
deterministic model is plausible but ad-hoc. 

3.1 Heating from Turbulent Current Sheets 

I begin with the view that the current sheets arise 
as the evolution of 2D, or quasi-2D MH D turbulence 



dZank and Matthaeusl 1 19921; ISpangler 1 1 19991) . In this case, 



the domain size A may be plausibly identified with the outer 
scale of the turbulence. The Joule heating in each current 
sheet £ is 



£ = r\ j V = rjj L c t c Z c 



(3) 



where V is the volume of a single sheet, given by V = 
L c t c Z c . The Joule heating from all the current sheets in the 
domain E is then given by 



E = N£ = N-qj 2 L c t c Z c 



(4) 



where N is the number of current sheets per domain. The 
mean volumetric heating rate in the domain e, which is taken 
to be the overall volumetric heating rate, is 



E 



Niqj 2 L c t c Z c 
A 2 A Z 



(5) 



where Vd — A 2 A 2 is the volume of a domain. Using the fact 
that the current per sheet is / = jL c t c , we have 



( Nr)Z c 
\A 2 A Z 



r- 



Lrtr 



(6) 



The question now arises as to how to relate the current 
in an individual current sheet, /, with the total current I b s 
within the Amperian Loop. This relation will depend on 
the model for the current sheets. For the remainder of this 
subsection, I will adopt the turbulence model in which there 
are, on average, equal numbers of positive and negative cur- 
rent sheets, and statistical fluctuations are responsible for 

lobs ? 0. 

Let Nt be the total number of current sheets within the 
Amperian Loop. We then identify the measured current I b s 
with the rms fluctuation in the total current contained within 
the Loop, 



lobs =Vn t I 

The total number Nt is given by 



A 

N T = ^N 



IS; 



las 



A 2 



N 



(7) 



(8) 



where A is the area of the Amperian Loop, I is the spacing be- 
tween the lines of sight, introduced in Section 1, and Si os is 
the effective line-of-sight extent of the coronal plasma. Equa- 
tion (8) is for the simplest case, in which the Amperian Loop 
is perpendicular to the large scale field. In the general case, 
a cosine of an orientation angle would be introduced in the 
numerator. This detail is ignored in the present discussion. 
Substitution of equations (7) and (8) into equation (6) yields 
the volumetric heating rate in terms of the measured total 
current I b s , 



V Z c 



I 2 

1 obs 

ISlns 



(9) 



As a final approximation, I assume that the extension of the 
domain and that of the current sheet along the large scale 
field direction are described by the same anisotropy index fi, 
A 2 = jLtA, Z c = fiL c . Use of these relations gives us the 
basic expression for the average volumetric heating rate due 
to turbulent current sheets in terms of the observed parameter 

lobs 



(1) 



1 



T 2 

1 obs 
ISlos 



(10) 



This expression factors itself neatly into three terms, each 
contained within brackets. The first is determined by the re- 
sistivity in the plasma and the domain properties. The second 
is determined by properties of the current sheets, specifically 
their thickness. The final term collects properties of the ob- 
servations, such as the inferred total current and the parame- 
ters of the lines of sight. 
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3.2 Heating from Deterministic Current Sheets 

In this subsection, I consider the possibility that the current 
sheets are not entirely random, and that there may be some 
preference for one sign of the current density, probably de- 
termined by the polarity of the large-scale coronal field. The 
total net current from the Sun must obviously be zero. How- 
ever, it is possible for a net current to exist in a limited region 
probed by a radio remote sensing measurement. We assume 
that the properties of the individual current sheets can be de- 
scribed as previously, so that the equations of Section 3.1 up 
to, and including equation (6), are valid. However, in the 
present case, there will be a different relationship between 
the total current I f, s and the current of an individual sheet, 
/. If there is a preference for current sheets of one sign of the 
current density, we can write 



bs 



(AT 4 



N-)I 



where N + is the number of sheets with positive current 
within the Amperian Loop, and A_ is the number of sheets 
with negative current density. The individual sheet current / 
is then taken as an absolute magnitude, with the sign of the 
current assumed in N + and A_ . If we introduce probabili- 
ties that the current densities will be positive or negative by 
N + = f+Nx, A_ = f-Nx, we have an expression for the 
current in a single current sheet, 



To. 



?,.s 



JV T (l-2/_) 



(12) 



The total current sheet number Nx is the same as that defined 
in equation (8). Substitution of equation (12) into (6), and 
algebraic manipulation gives the volumetric heating rate in 
the case of "deterministic" current sheets 



~N 



1 



*c(l-2/-) 



I 2 

1 obs 

i 2 sL 



(13) 



where the expression has again been factored into terms 
which contain, respectively, characteristics of the plasma, the 
current sheets, and the observations. 

3.3 Comparison of the Expressions for the Heating Rate 

The expressions for the volumetric heating rate in the two 
models of the current sheets, equations (10) and (13) respec- 
tively, appear quite different in form, and it is natural to ask 
which is the larger for realistic input parameters. In other 
words, given a measurement of I b s , would greater Joule 
heating result if the current were distributed in a random set 
of turbulent current sheets as described in Section 3.1, or in 
a set of sheets with predominantly one sign of the current 
density, as discussed in Section 3.2. 

Let the heating rate expression for a turbulent set of current 
sheets as given in equation (10) be noted by ex, and that due 
to a systematic set of sheets with preferentially one sign of 



the current density (equation (13))as eg. Equations (10) and 
(13) can be easily manipulated into the following form 



(_s_ 



1 



iV(l-2/_) s 



/ A 2 



\lSlo 



(14) 



If one assumes that the first term in square brackets on the 
right hand side of this equation is of order unity, then the 
relative heating rate depends on the ratio of the domain area 
to that of the Amperian loop. The precise value of this ratio 
depends on the circumstances of the observations, as well as 
the value of A. An estimate of its value in the case of the 



observations of Spa ngler I (12007b is given at the end of the 
next section. In what follows, I discuss the case of turbulent 
current sheets, then briefly note that the conclusions would 
not be significantly different for the deterministic case. 



(11) 3.4 Estimate of the Turbulent Heating Rate 



Equation (10) is now used to estimate the coronal heating rate 
from turbulent current sheets. The variables in the last term 
(lobs, I, Sios) are observational parameters and are known. 
The calculation will be carried out for the conditions charac- 
teristic of the large ARM event of August 16, 2003, in which 
the inferred current was 2.5 x 10 9 Amps. The line of sight 
to the radio source had a minimum heliocentric distance of 
6.7i?0. Coronal plasma properties characteristic of this dis- 
tance will be used in the calculation below. A similar a nalysis 
at other times in the observations of I Spangler" ( 20071) . when 
there were only upper limits to the current, would obviously 
yield lower values for the Joule heating rate. The calculation 
also requires estimates of 77, A, and t c . 

3.4.1 Resistivity 

For the resistivity 77, the Spitzer resistivity is used, which is 
based on Coulomb collisions of current-carrying electrons 
with ions and other electrons. It is certain to be a drastic 
underestimate, in that the true resistivity is almost certainly 
determined by collisionless processes. However, the Spitzer 
resistivity can be derived from fundamental principles, which 
is not true of other estimates, and it can serve as a lower limit 
to the true resistivity. An informal discussion of the possible 
role of collisionless processes in determining the resistivity 
is given in Section 4 below. 

The S pitzer resistivity is the reciprocal o f the conductivity 
given bv lGurnett and Bhattachariee I (|2005) 



32^e 2 (2fc B r e ) 3 / 2 
z 2 lnA 



(15) 



In this equation, and equation (16) below, A stands for the 
Coulomb logarithm rather than the domain size as used oth- 
erwise. The electron temperature is T e . All other terms in 
equation (15) have been defined, or are obvious fundamental 
physical constants. 



6 



Steven R. Spangler: Currents in Solar Corona 



Equation (15) 
resistivity in a 



can be used 
"suitable for 



to write the Spitzer 
observers" form as 



dGurnett and Bhattachariee L |2005 ) 
In A 



r]s = 5.2 x 10 



-5 



Ohm-m 



(16) 



where the thermal energy ksT e is now given in elec- 
tron volts. For approximate coronal conditions I choose 
a value for the Coulomb logarithm of A = 25 



(Kral l and Trivelpiece L [1973). With an assumed coronal 



temperature of 2 x 10 6 K, appropriate for closed-field re- 
gions (electron thermal energy ksT = 172 eV in equation 
(16)), the resistivity is i~is — 5.74 x 10~ 7 Ohm-m, or about 
35 times the resistivity of silver. 

3.4.2 Domain Size 

I will take the domain size A to be the outer scale of the 
turbulence in the relevant part of the corona. There are 
two estimates in the literature for this outer scale. The 
first is mean spacing between flux tubes whi ch expand into 



the co rona. This estimate was intro duced b y Hollweg et al 



1982), an d subsequently used by Mancuso an d Spangler 



1999) and Cranmer and van Ballegooiien (2005). The for- 



mula used bv lMancuso and Spangler Id 19991) is 



1.37 x 10 7 
A = — meters 

B(G) 



(17) 



where B(G) is the magnetic field strength in Gauss. For the 
magnetic field in th e corona, we use the recent estimate of 
Inglebv et al 1 d2007l) which was obtained f rom Faraday rota - 
tion measurements very similar to those of I Spangler! d2007t) . 
They found that the magnetic field could be represented by 
an inverse square dependence on the heliocentric distance, 
with a normalizing value of ~ 0.050 G at r = 5Rq. At 
a heliocentric distance of 6.7 Rq, the estimated magnetic 
field is 2.78 x 10 -2 G, and the corresponding value of A 
is 8.2 x 10 7 m. In a more recent theoret ical study of c orona l 
heating and solar wind acceleration, ICranmer et al ] d2007l) 
argue for a smaller value of the domain size (their parameter 
L± which serves as the outer scale of the turbulence), which 
physically corresponds to the diameter of the photospheric 
flux tubes rather tha n their separation. A dopting the estimate 
for A = L± from ICranmer et al I (120071) would reduce our 
value of A by about a factor of 4 from the estimate of equa- 
tion (17). 

A second estimate of the outer scale of coronal turbulence 
comes from power spe ctra of fluctuating Doppler shifts of a 



space craft transmitter ( Wohl muth et al L 2001; E fimov et al 
2004). These estimates, which result from measurements 



rather than plausible theoretical arguments, give outer scales 
from a few tenths of a solar radius to a solar radius or more 
at heliocentric distances of 5 — 10-R^. The values reported 



bv I Wohlmuth etall d200lh and lEfimov et all d2004 are sev- 
eral times larger than that given by equation (17). It is ob- 
vious that the factor of 4 smaller value for A advocated 



by ICranmer et al I Q2007) is in more serious disagree ment 
with the observatio nal value of I Wohlmuth et al I (120011) and 
Efimo v et al I d2004 . A resolution of this matter would war- 
rant a paper in its own right, but for the present work we 
use equation (17). As may be seen from the heating rate ex- 
pression in equation (10), lower values of the domain size A 
generate higher values of the heating rate e. 

3.4.3 Current Sheet Thickness 

For the current sheet thickness t c , I choose the ion inertial 
length t c = where Va is the Alfven speed, and f2j is 
the proton ion cyclotron frequency. This would seem to be 
both plausible and a good lower limit to what the current 
sheet thickness can be. Once again, equation (10) shows 
that use of a minimum plausible value for t c leads to an up- 
per limit to the heating rate e. To calculate the ion inertial 
length, the plasm a density profile given by equation (6) of 
I Spangler I (12007) (based on radio propagation measurements 



of the corona) and the magnetic field model of llngleby et al 



(2007) are used. These yield the following formula for the 
estimated current sheet thickness 



Rc 



/-,. = 1 .0 x 10 J I | meters 



1.2 



(18) 



where Rq is the heliocentric distance in units of a solar ra- 
dius. For R = 6.7i? , t c = 1.4 x 10 3 m. 

3.4.4 Observational Parameters 

The observed parameters in equation (10) are contained in 
the term in the third set of brackets. The observed current 
lobs — 2.5 x 10 9 A mps and the sepa ration of the lines of 
sight I is 33,000 km dSpangler L 120071) . There remains the 
value for the effective thickness of the plasma alon g the line 
of sight. I use the expression from lSpanglerl (2002) 



Si, 



-RoR^ 



(19) 



Use of the above parameters with an impact parameter 
i?o = 6.7 gives a heating rate e = 1.27x 10~ 16 Watts/m 3 . To 
determine the signific ance of this number, I compare i t to the- 
ore tical estimates of ICranmer and van Ballegooijen I d2005l) 
and lCranmer et al ( 20071) . These papers utilize cgs units, and 
report heating rates in power per unit mass. The volumetric 
heating rate given above is then 1.27 x 10~ 15 ergs/cm 3 /sec, 
and is converted to a heating rate/unit mass q 



q = e/p = 3.0 x 10 4 ergs/sec/gm 



(20) 



where I have again used the power law density model of 
equation (6) of I Spangler I (120071) in obtaining the mass den- 
si ty at r = 6.7i?p . 



Cranmer et al I d2007l) calculate heating rates as a function 
of heliocentric distance. Since their calculations are self- 
consistent, their heating rates may be considered to be those 
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which are required by the observed heating of the corona 
and acceleration of the solar wind. Examination of Fig- 



ure 8 of ICranmer et al I (2007) shows that the heating rate 
per unit mass at a heliocentric distance of ~ 7Rq is in the 
range of 10 10 — 10 11 ergs-sec _1 -gm _1 , depending on the as- 
sumed amplitude of the photospheric velocity fluctuations. 
We therefore conclude that the heating rate given by equa- 
tion (10) is lower than values which are required to account 
for coronal heating by at least a factor of 3 x 10 5 , if the input 
parameters used here are valid. If the outer scale to the turbu- 
lence is a factor of ~ 4 le ss than the value given by equation 
(17), as recommended by lCranmer et al ( 20071) . the ratio of 
mass heating rates would be about 10 5 . In either case, this 
huge mismatch means that exercises with fine tuning the pa- 
rameters in the model would be a fool's errand. It should be 

A 2 

noted that the ratio ^ — which appears in equation (14) is 
of the order of unity, within a factor of several either larger 
or smaller depending on the assumed outer scale of the tur- 
bulence. The conclusion on the magnitude of Joule heating 
would not be changed by adopting the non-turbulent current 
sheet model of Section 3.2. 

There are two possible conclusions to be drawn from the 
calculations of this section. 

1. In view of the large disparity between the calculated 
Joule heating rate and that which is required for a signif- 
icant contribution to the thermodynamics of the corona, 
the currents which may have been observed are irrele- 
vant for coronal heating. This argument would seem to 
be strengthened by the fact that I used the largest de- 
tected value of lobs from the two days of observation. 
Other intervals would have provided smaller values for 
lobs or upper limits thereto, yielding smaller values of 



2. A more likely explanation, in my opinion, is that these 
current systems do play an important role in coronal 
heating, but that role is underestimated in the calcu- 
lations presented here, because they are based on the 
Spitzer resistivity. According to this viewpoint, the 
analysis of Sections 3.1 and 3.2 is valid, but a cor- 
rect calculation would require an appropriate, and much 
larger value of the resistivity. This is clearly speculation 
until it can be demonstrated that a much larger resistiv- 
ity (by orders of magnitude) characterizes the coronal 
plasma at 5Rq < r < IORq. 



4 Possible Enhancement of Resistivity in Current 
Sheets 

In this section, I consider the second of the possibilities listed 
immediately above, i.e. that the resistivity could be suf- 
ficiently enhanced in these coronal current sheets to make 
Joule heating a thermodynamically important process. An 
obvious way for this to happen is a plasma instability that 



produces high levels of fluctuating electric or magnetic fields, 
which scatter the current-carrying electrons and enhance the 
resistivity. To assess this possibility, we need to examine 
the magnitude of the electron drift speed within the current 
sheets. 

Equations (7) and (8) give the relationship between the ob- 
served current in the Amperian loop, I bs, and the current in a 
single sheet, I. Using these equations and the identity imme- 
diately before equation (6), we have for the current density 
in a single sheet 



J 



Mo 



t c L c ^f{lSios)^fN 



and for the electron drift speed 



Al 



vd 



en ent c L c ^f{lSios)\fN 



(21) 



(22) 



where e is the fundamental electric charge. To simplify equa- 
tion (22) I adopt a set of plausible assumptions. I assume that 
the width of a current sheet will be some fraction of the do- 
main size, L c = (3 A with probably having a value between 
0.1 and 0.5. From Section 3, we already have estimates of 
other parameters (e.g. t c , L c , n) in equation (21) at the fidu- 
cial heliocentric distance of 6.7Rq. This yields the following 
estimate for the electron drift speed, 



2.17 x 10 6 , 
vd = 1= — m/sec 



(23) 



For this expression to be meaningful in the context of 
plasma instabilities, we need to compare it with a charac- 
teristic plasma speed. An obvious choice is the electron ther- 
mal speed vg — J ( Nicholson , 1983 ). I use a value of 

T e = 2 x 10 6 K, which is characteristic of closed-magnetic- 
field regions in the corona. Open field regions would have 
a lower temperature and lower thermal speed. We then have 
for the drift speed to thermal speed ratio 



V0 



0.39 



(3VN 



(24) 



As mentioned in the definition of (3 immediately above, 
and the discussion of N in Section 3, (3 is a number which is 
probably less than unity, but not by a large factor, and N 
is an integer which is probably larger than unity, but not 
much greater. Their product should therefore be of order 
unity. This calculation then suggests the plausibility of elec- 
tron drift speeds of order the thermal speed in these current 
sheets. If the drift speed is of the order of the electron ther- 
mal speed, it is also of order or larger than the ion acoustic 
speed. 

It is necessary to stress that this calculation has contained 
products of several parameters (such as A, L c , etc) which are 
imperfectly known, so the net result presented here is sim- 
ilarly uncertain. However, the conclusion of this section is 
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that a current-driven instability, which would lead to high 
levels of fluctuating electric and magnetic fields, is a possi- 
bility. 

This calculation has also been carried out for the case of 
"deterministic" current sheets discussed in Section 3.2. The 
details of the calculation are not presented here, but the fi- 
nal result is that the drift-to-thermal speed ratio is somewhat 
smaller (a factor of about 0.17), but not enough to alter the 
qualitative conclusion stated above. 

The possible existence of substantial electron drift speeds, 
comparable to the electron thermal speed, raises the possibil- 
ity of an interesting obs ervational diagnostic of such current 
sheets. I Spangler I d 1998b pointed out that an electron distribu- 
tion carrying a current will have its distribution function dis- 
torted in the direction of current flow, and accordingly have a 
more populated tail than a distribution which carries no cur- 
rent. This additional tail component makes the plasma more 
effective at collision excitation of ions to excited states whose 
energy abo ve the ground stat e is a few times the electron ther- 
mal energy. ISpangler I dl998D suggested that turbulent current 
sheets might reveal themselves via enhanced emission line 
glow as the excited ions radiatively de-excite. An important 
parameter determining the intensity of the line radiation is 



A, 



2k B T 



(25) 



This parameter is approximately equal to the square of the 
drift speed to thermal speed ratio. When Ad becomes of 
the order of a few tenths, the line emission can be substan- 
tially ^enhanced rel^veto_flie_£urrenfTree value (see Figure 
9 of ISpanglerLfl998l) . ISpanglerl dl998l) found that A D « 1 
for turbulence in the interstellar medium, so turbulent en- 
hancement of emission line radiation probably does not occur 
there. However, the results presented in this section suggest 
that this mechanism is much more likely to occur in the solar 
corona. 

4. 1 The consequences of current-driven instabilities 

A complete discussion of the consequences of a current- 
driven plasma instability for the resistivity within coronal 
current sheets is beyond the scope of the present paper. I 
will only briefly refer to some results in the literature which 
indicate the effect may be significant. 

The issue of in stabilities d ue to high electron drift speeds 
was discussed i n Spangler (1998), who cited results from 



Drummond and Rosenbluth (1962). The remarks made 
there are still relevant to the present discussion. The sum- 



mary of the work of Drummo nd and Rosenbluth I (119621) pre- 



sented in I Spangler I (11998b is that an electron drift speed 
greater than ~ 0.1 2vg could be sufficient for excitation of 
obliquely-propagating electrostatic ion cyclotron waves. 
The role of a current-dri ven instabi l ity in e nhanced re 



sistivity was discussed by IChittenden I (119951) . I Chittenden 



laboratory Z-pinches which showed these structures to be 
larger than expected on the basis of theory with a Spitzer 
conductivity. The unexpectedly large size of Z pinches 
suggested that enh anced transport coefficients were present. 
Chittenden I (119951) found that the electron drift speed ex- 
ceeded the ion acoustic speed in the outer edges o f the Z 
pinch. His results (see Figure 3 of IChittenden I ((.1995)) 
showed that the effective resistivity due to lower hybrid 
waves could exceed the Spitzer resistivity by 3 - 4 orders of 
magnitude. This enhancement is approaching that needed for 
thermodynamic relevance of coronal currents, as discussed in 
Section 3.4. 

Enhanced resistivity has also been hypothesized to play an 
important role in magnetic reconnection, allowing reconnec- 
tion to proceed at a faster ra te and produce heating of the 
plasma. iKulsrud et al I (12005) discussed experiments show- 
ing the presence of magnetic fluctuations wit hin the recon- 



nectio n current sheet on the MRX experiment. IKulsrud et al 



(2005) identify these fluctuations as obliquely -propagating 
waves arising due to a cross-field curr ent which has a drif t 
speed exceeding the Alfven speed. IKulsrud et al ] d2005l) 
found an enhancement of the resistivity, estimated from the 
wave force on the electrons, which exceeds the Spitzer resis- 
tivity by a factor of several. The relatively modest enhance- 
ment of the resistivity relative to t hat required for co ronal 
relevance, or the results discussed bv lChittenden (1995), can 
be attributed to the nearly collisional dynamics of the MRX 
experiment. 

The discussion in Kulsrud et al J 20051) has recently been 



super ceded by the results of Wa ng. Kulsrud. and Ji I (2007 



2008), who now favor perpendicularly-propagating, unsta- 
ble waves which nonlinearly couple to magnetosonic waves. 
It is the magnetosonic waves which determine the resistiv- 
ity. This anomalou s resistivity is higher than estimated in 
Kulsrud etalld2005l) . 

The above-cited studies are not intended to correspond in 
detail to the case of Joule heating of the solar corona and the 
highly enhanced resistivity tha t would be required the r e. The 
inve stigations of IChittenden I (fl995h. IKulsrud et al I d2005l) 
and IWang. Kulsrud. and Ji I d2007i 120081) are of importance 
in showing that current filaments with drift speeds compara- 
ble to or exceeding the ion acoustic speed can produce wave 
and turbulence fields that substantially enhance the resistiv- 
ity. Current sheets or filaments with drift speeds approaching 
the electron thermal speed would be even more subject to in- 
stability, and to a wider range of unstable modes. 

A final point to be considered in this section is whether the 
resistivity in the coronal current sheets could plausibly reach 
levels necessary for important Joule heating. The arguments 
in the previous paragraphs have shown that current-driven 
instabilities could quite plausibly be present, but could they 
enhance the resistivity by several orders of magnitude? A 
very general expression for the resistivity is 



d 19951) developed a fluid theory to explain observations of 



(26) 
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where v is a collision frequency of some sort. If collisionless 
scattering of electrons determines the resistivity, v could be 
as large as the frequency of a high frequency plasma mode. 
In what follows, I choose the lower hybrid frequency 



Z7T 



(27) 



In this identity, Hi and il e are, respectively, the ion and elec- 
tron cyclotron frequencies. 

I use the lower hybrid frequency as a proxy for a plasma 
mode frequency which is above the ion cyclotron frequency, 
and do not claim that the unstable waves in coronal current 
sheets are nec essarily lower hyb rid waves. In support of 



this approach. IChittenden ( 1995b presents estimates of the 



anomalous collision frequency which are of order ujlh, with 
a multiplicative constant dependent on t he drift speed to io n 
acoustic speed ratio (see equation (2) of Chittenden , 1995) 



Usin g the coronal magnetic field model of llnglebv et al 



d2007h . we estimate a lower hybrid frequency vlh = 1-82 
kHz at a heliocentric distance r = 6.7 Rq. Substitution of 
this collision frequency into equation (25) (again using the 
same estimate of the electron density n used in equation 
(22)) gives an anomalous resistivity of 6.27 Ohm-m. This 
exceeds the Spitzer resistivity calculated following equation 
(16) by approximately 7 orders of magnitude. This en- 
hancement in the resistivity is comparable to, and in fact 
exceeds, the factor by which our calculated heating rate 
must be increased in order to be relevant for the thermody- 
namics of the corona dCranmer and van Ballegooijen 



Cranmer et al L 120071) . This brief calculation then suggests 



that self-enhanced resistivity in coronal current sheets could 
lead to thermodynamically -relevant levels of Joule heating in 
the corona. 



in the coronal plasma, and could show a preference for 
one sign of the current density. These derivations are 
given in Sections 3.1 and 3.2, and provide the formulas 
for the volumetric heating rates given in equations (10) 
and (13). 

3. Use of th e se for mulas, with observational data from 
Spangler] (2007) and plausible independent coronal 
data, yield an estimated heating rate of 1.3 x 10 -16 
Watts/m 3 (1.3 x 10~ 15 ergs/cm 3 /sec in cgs units). The 
corresponding heating rate per unit mass is 3.0 x 10 4 
ergs/gm/sec. This appears to be smaller than the level 
necessary to be significant for coronal heating by at least 
a factor of 3 x 10 5 . 

4. The conclusion to be drawn from point (3) is that ei- 
ther these currents are irrelevant for coronal heating, or 
that the true resistivity in the corona exceeds the Spitzer 
value by several orders of magnitude. Resolution of this 
matter obviously lies in a better understanding of the re- 
sistivity in a collisionless plasma. 

5. The same model used to estimate the volumetric and 
mass heating rates is also used to estimate the electron 
drift speed in the current sheets. This drift speed could 
be comparable to the electron drift speed, and in excess 
of the ion acoustic speed. Accordingly, current-driven 
instabilities might be present in these sheets, and the 
waves driven unstable by these currents might enhance 
the resistivity to significant levels. This contention is 
supported by works in the literature which have shown 
enhancement of resistivity by current-driven instabili- 
ties. 



5 Conclusions 



1 . Radioastronomical observations reported by I Spangler 



(2007) are consistent with coronal currents flowing 
through Amperian loops defined by adjacent lines of 
sight to different components of a radio source. There 
are estimates of currents of 2.5 x 10 9 and 2.3 x 10 8 
Amperes, respectively, on two days. Another interval of 
high quality data on one of the days yielded an upper 
limit to the differential Faraday rotation, and a corre- 
sponding upper limit to the current of 8 x 10 s Amperes. 
These data are used as input for a calculation of Joule 
heating of the solar corona. 

2. Two models are developed to calculate the Joule heating 
associated with the observed currents. In both models, 
the current is envisioned as being in thin, intense cur- 
rent sheets stretched out along the large-scale coronal 
magnetic field. The first views the sheets as arising in 
the evolution of quasi-2D magnetohydrodynamic turbu- 
lence. The other assumes that current sheets will arise 
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